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ABSTRACT 


Experiments  with  a  20mm  M 56  shell,  whose  cylindrical  cavity, 
fineness  ratio  2.7,  was  filled  with  different  amounts  of  various 
liquids,  demonstrated  the  existence  of  dynamic  instability  for  a 
broad  range  of  loading  conditions.  It  was  found  that  for  water¬ 
like  liquids,  the  nutational  component  of  yaw  is  divergent  when  40 
percent  to  90  percent  of  the  cavity  is  filled.  For  mercury,  the 
yaw  is  divergent  from  about  25  percent  to  100  percent.  The  severity 
of  instability  or  the  rate  of  divergence  of  yaw  increases  with  the 
increase  in  the  specific  gravity  of  the  liquid. 


Stewartson’s  theory  applies  only  to  the  steady  state  when  the 
liquid  is  spinning  with  the  full  spin  of  the  shell.  Under  this 
condition  the  theory  predicts,  for  a  cylindrical  cavity,  at  what 
loading  conditions  (per  cent  of  fill)  the  liquid  filled  shell  system 
will  be  unstable.  Predicted  instabilities  are  usually  confined  to 
relatively  narrow  zones  of  loading  conditions.  The  observed,  very 
broad  zones  of  instabilities,  are  attributed  to  the  transient  state 
of  liquid  rotation.  The  time  required  to  reach  the  steady  state  is 

ftti  -i  rol  *\r  1  nnnr  oc*nor»*5  al  lir  1  imvirlc  t.H  +Vi  -\r^  onnoi  -H  ac  1  coo  4~Vi  nn 
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say,  10  centi-stokes  (c.s.).  During  this  time,  if  the  shell  is 
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unstable,  the  yaw  may  reach  a  prohibitively  large  level  rendering 
such  shell  useless  in  practice.  Baffles  in  the  cavity  effectively 
suppress  the  transition  phase. 
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1.  INTRODUCTION 


It  is  a  well-established  fact  that  some  spin  stabilized  shell 
if  filled  with  a  liquid  become  dynamically  unstable.  The  problem  is 
well  illustrated  by  Plate  1.  The  cylindrical  cavity  of  this  particular 
20mm  shell  was  filled  with  mercury.  The  shell  is  gyroscopically  stable 
with  a  stability  factor  of  about  2.5-  It  began  its  flight  with  very 
small  yaw  but  at  280  feet  from  the  muzzle  the  yaw  has  increased  to  about 
80°.  Such  exceptionally  rapid  divergence  of  yaw  is  usually  associated 
with  heavy  liquids.  Lighter  liquids  cause  milder  instabilities.  Never¬ 
theless,  even  milder  instability  persisting  for  a  sufficiently  long 
time  may  render  the  shell  useless  in  practice. 

Although  the  problem  of  stability  of  liquid  filled  shell  has  been 
under  investigation  for  a  long  time,  unfortunately,  to  date,  no  infallible 
general  principle  applicable  to  all  cases  has  emerged. 

Largely,  through  experience  supplemented,  however,  by  theoretical 
analyses,  the  following  variables  appear  to  influence  the  dynamics  of 
liquid  filled  shell: 

1.  Air  space  in  cavity  or  per  cent  of  filled  volume. 

2.  Geometry  of  the  cavity. 

3.  Specific  gravity  of  the  liquid. 

4.  Liquid  viscosity. 

5.  Spin  level  of  the  shell. 

6.  Shell  velocity. 

7*  Quadrant  elevation. 

It  is  not  yet  possible  to  specify  the  relative  importance  of  these 
variables  in  general.  However,  for  a  cylindrical  cavity,  the  air  space 
and  the  geometry  of  the  cavity  are  probably  the  most  important,  followed, 
perhaps,  by  the  specific  gravity  of  the  liquid.  The  viscosity  has  at  least 
two  effects.  First,  higher  viscosity  expedites  the  attainment  of  full  spin 
and,  hence,  shortens  the  duration  of  the  transitional  regime  which,  as  we 
shall  see  later,  is  very  important.  Second,  higher  viscosity  probably  in¬ 
hibits,  within  the  disturbed  fluid,  the  development  of  oscillations  of 
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large  amplitude  and  causes  these  disturbances  to  die  out  faster.  The  spin 
level  of  the  liquid,  rotating  as  a  rigid  body,  controls  the  frequencies  of 
its  free  oscillations.  Since  instability  of  the  shell  occurs  only  whenever 
one  of  the  liquid  frequency  is  in  resonance  within  the  nutations!  frequency 
of  the  shell,  the  spin  of  the  shell  is  important.  It  can  be  shown  (see  W.  £. 
Scott's  forthcoming  BRL  report)  that  the  frequencies  of  free  oscillations 
of  the  liquid  are  given  by 
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where  SI  is  the  spin  of  the  liquid,  S2  is  that  of  the  shell,  and  f  is  a 
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certain  function  which  defines  the  double  infinity  of  frequencies  of  the 

liquid.  The  nutational  frequency  of  the  shell  is  given  by 
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where  A,  B  are  the  polar  and  transverse  moments  of  inertia  of  the  shell  and 
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where  s  is  the  gyroscopic  stability  factor.  This  factor  varies  as  ■  s  • 

0 

m 

a 

The  C  ,  the  overturning  moment  coefficient,  depends,  among  other  things, 


m 

a 

on  Mach  number  or  velocity.  Thus,  the  spin  and  the  velocity  change  a  and, 
hence,  the  nutational  frequency.  Theoretically,  therefore,  it  is  possible, 

with  the  liquid  spinning  with  the  full  spin  or  the  shell  —  =  1,  to  alter 
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spin  or  velocity  (through  change  in  C  )  or  both.  The  effect  of  quadrant 
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elevation  for  higii  angle  fire  may  well  be  ttasociauea  wii^n  cnangmg  re±a- 

^i  i 

tionship  between  the  — — 1 —  of  the  liquid  and  the  nutational  frequency  of  the 
shell.  e 


In  the  history  of  studies  of  this  problem,  two  milestones  should 
be  briefly  mentioned.  In  19^0,  E.  A.  Milne^"  made  a  thorough  theoretical 
analysis  of  this  problem  on  the  basis  of  a  mathematical  model,  developed 
earlier  by  others,  of  the  stability  of  a  shell  containing  a  spheroidal 
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cavity  completely  filled  with  inviscid  fluid,  rotating  with  the  full  spin  of 
the  shell.  The  stability  of  such  a  system  depends  on  the  geometry  of  the 
cavity,  the  major  and  minor  axes  of  the  spheroidal  cavity,  the  mass  of  the 
liquid,  and  the  inertial  properties  of  the  shell.  He  showed  that,  for  this 
model,  the  region  of  instability  was  defined  by  the  following  inequalities : 

x>y,  x-y<cl<x  +  y 

where 


where  K  =  fineness  ratio  of  the  cavity,  ^ 

M  =  mass  of  the  liquid 

B,  A  =  transverse  and  axial  moments  of  inertia  of  the  shell. 

When  all  the  firings  from  1926  onward  were  analyzed  and  each  shell 
was  marked,  according  to  its  behavior,  as  either  S  =  stable  or  U  =  unstable, 
and  plotted  in  the  (x,  y)  plane,  it  was  found  that  a  fairly  well-defined  curve 
could  be  drawn,  in  this  plane,  separating  the  S  from  the  U  points.  This  curve, 
v>ATj0Tror  D-nnoarorl  +0  Vw=>  nnlv  t.Ancr^nt.  to  t.h p  theoretical  instabilitv  boundary 

uwnu  y  wv-*-  v\a  vv#  ^  ^  v**-*-j  ^  w.  — - - - -  %j  - %/ 

x  -  y  =  1  near  y  =  0  and  was  purely  empirical  elsewhere.  Nevertheless,  it 
proved  to  be  very  useful  in  practice. 


Milne  further  explored  the  effect  of  viscosity  and  of  air  space  and 
suggested  how  these  parameters  would  influence  the  stability.  However, 
such  effects  could  not  be  adduced  quantitatively  from  the  available  theory. 

2  3 

The  next  high  light  in  this  story  is  the  work  of  K.  Stewartson  ’ 

He  extended  the  analysis  to  a  cylindrical  cavity  partially  filled  with  an 
inviscid  fluid  rotating  as  a  rigid  body  with  the  full  spin  of  the  shell. 

This  model,  of  course,  is  much  more  realistic  than  the  completely  filled 
spheroidal  cavity  but,  also,  is  considerably  more  complicated.  Stewartson 
showed  that  instability  results  only  when  the  poles  of  a  certain  very  complicated 
nondimens ional  function,  defined  entirely  by  the  dimensions  of  the  cavity  and  the 
fluid  mass,  occur  close  to  the  nutational  frequency  of  the  shell.  In  certain 

■P4 n  4--v*d  Q+  TDr\  -n-h  +-  -?  0V1  PT*mHnrr  rir*oirnr}  -hVr1«  thpArv  r  *=>  rl  « 
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appeared  to  be  fully  confirmed.  For  different  fineness  ratios  of  the 
cylindrical  cavities  and  different  air  spaces,  all  shell  predicted  to  be 
unstable  gave  short  ranges  and  all  shell  predicted  to  be  stable  gave  full 
range.  This  problem,  therefore,  appeared  to  have  been  solved. 

Recently,  however,  at  the  Aberdeen  Proving  Ground,  in  the  course  of 
development  of  a  certain  weapon  system,  we  have  encountered  a  severe 
case  of  dynamic  instability  of  liquid  filled  shell.  For  the  geometry  of 
its  cylindrical  cavity  and  for  various  air  spaces  which  have  been  fired, 
it  should  have  been  stable  according  to  Stewartson's  analysis,  but  it 
wasn't.  It  was  violently  unstable  under  all  conditions  of  test  firings. 

Further  search  of  the  literature  proved  not  very  helpful.  A  very 
extensive  amount  of  both  theoretical  and  experimental  work  on  sloshing 
of  fluids  in  tanks,  in  connection  with  the  development  of  liquid  fueled 
rockets,  has  been  done  in  recent  years.  There  is  a  great  variety  of 
theoretical  treatments  of  related  problems  but  mostly,  however,  with 
nonspinning  fluids.  A  fairly  extensive  bibliography  of  papers  in  this 
field  is  to  be  found  in  reference  (4). 

Therefore,  it  appeared  desirable  to  look  into  this  problem  further. 
Since  detailed  experimental  data  on  the  dynamics  of  liquid  filled  shell 
are  practically  nonexistent,  it  appeared  particularly  desirable  to  obtain 
such  information.  Availability  of  the  spark  ranges  at  BRL  made  the 
acquisition  of  such  data  feasible. 

On  the  theoretical  side,  we  set  out  to  examine,  within  the  frame¬ 
work  of  modem  ballistic  theory,  four  fundamental  cases.  In  the  order 
of  increasing  complexity,  these  are: 

1.  Nonspinning  liquid 

a.  Completely  filled  cylindrical  cavity 

b.  Partially  filled  cylindrical  cavity 

2.  Spinning  liquid 

a.  Completely  filled  cylindrical  cavity 

b.  Partially  filled  cylindrical  cavity 
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Problems  la  and  lb  have  been  solved  and  the  reports  published  by- 
Scott,  references  (5)  and  (6).  While  some  of  the  conclusions  are  not 
new,  Scott's  independent  formulation  and  analysis  of  these  problems 
show  much  more  clearly  than  we  were  able  to  find  in  the  available 
literature,  the  effect  of  the  liquid  on  the  yawing  motion  of  the  shell. 

His  work  is  particularly  valuable  for  the  comparison  of  the  theoretical 
predictions  with  the  observed  behavior  of  the  shell  in  flight.  The 
available  treatments  of  these  problems  require  a  considerable  amount  of 
adaptation,  and  are  also  either  too  restrictive  or  too  general  to  be 
directly  useful. 

Scott's  analysis  shows  that  for  a  completely  filled  cavity,  the  only 
effect  the  liquid  produces  on  the  dynamics  of  the  shell  is  to  modify  its 
inertial  properties.  For  a  partially  filled  cavity  (nonspinning  liquid), 
the  liquid,  in  addition  to  altering  the  inertial  properties  of  the  empty 
shell,  affects  its  yawing  frequencies.  There  is  no  contribution  by  the 
liquid  to  the  damping  of  the  system.  Since  nonspinning  liquid  can,  in 
reality,  occur  only  in  nonspinning  shell,  the  conclusion  of  this  analysis 
is  that  the  liquid  filler  has  no  detrimental  effect  on  the  dynamics  of  fin 
stabilized  shell.  This  conclusion  has  been  verified  by  a  limited  number 
of  experiments^. 

Problems  2a  and  2b  have  been  solved  by  Stewartson.  Scott,  however, 
has  removed  Stewartson* s  assumption  of  small  liquid  mass.  For  heavier 
liquids,  this  modifies  the  nutational  frequency  of  the  system  and,  hence, 
changes  Stewartson’s  stability  regions.  Moreover,  by  carrying,  in  the 
analysis,  the  ratio  of  the  spin  of  the  liquid  to  that  of  the  shell  (  in 
Stewartson 's  analysis,  this  ratio  is  unity)  and  choosing  a  different  boundary 
condition  which  is  to  be  made  homogeneous,  Scott's  equations  are  more 
physically  understandable.  They  show  clearly  how  the  spin  of  the  liquid 
modifies  the  nonspinning  case.  The  analysis  of  problem  2a  is  being  published 
and  that  of  2b  will  be  published  shortly.  The  present  report  deals  only 
with  the  experimental  results. 
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2.  EXPERIMENTS 


The  test  vehicle  selected  was  the  20am  M56  (T  282)  spin  stabilized 
shell.  It  was  cheap,  abundant,  and  of  convenient  size  for  the  work  in  the 
spark  range.  Its  schematic  drawing  is  given  in  Figure  1.  The  normal 
cavity  of  this  shell  has  thicker  walls  under  the  rotating  band  because 
this  shell  is  designed  to  be  fired  at  high  velocity.  The  cavity,  therefore, 
was  nx>dified  to  make  it  cylindrical  with  a  fineness  ratio  2.68.  Most  of 
the  firings  were  done  from  a  gun  with  a  twist  of  rifling  of  one  turn  in 
25  calibers  of  travel  with  a  muzzle  velocity  of  about  27OO  fps.  The  shell 
were  fired  thru  a  spark  range.  This  range  is  fully  described  in  reference 
(7).  Briefly,  it  contains  45  spark  stations  over  a  base  line  of  285  feet. 
The  first  observing  station  is  15  feet  from  the  muzzle.  The  observed 
trajectory,  therefore,  is  about  4400  calibers  long  over  which  the  20mm 
shell  executes  about  15  yawing  periods.  The  accuracy  with  which  the  shell 
position  and  attitude  can  be  measured  on  photographs  is  0.01  inches,  and 
for  these  relatively  short  shell,  5  minutes  of  arc  respectively.  If  dis¬ 
similar  pins  are  inserted  in  the  base  of  the  shell,  their  orientation 
at  each  station  can  be  determined  to  an  accuracy  of  about  2  degrees. 

These  pins  are  used  to  measure  the  axial  spin  of  the  shell.  Thus  the 
range  data  provide  a  fairly  accurate  history  of  the  yaw,  the  axial  spin 
and  the  velocity  of  the  shell  over  the  initial  4400  calibers  of  travel. 

Three  pieces  of  data,  therefore,  are  directly  obtainable  from  the  photo¬ 
graphic  plates:  (a)  position  of  the  shell  at  a  given  time,  (b)  the 
attitude  of  the  shell  axis,  or  the  yaw,  and  (c)  the  orientation  angle 
of  the  base  pins  vector. 

Analysis  of  these  data  is  handled  in  the  usual  manner  which  is  fully 

described  in  reference  (8).  The  analysis  of  yaw,  for  example,  is  done 

by  fitting  an  epicycle  to  the  observed  components  of  yaw.  With  rigid 
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filler,  an  epicycle  fits  the  yawing  motion  of  this  shell  very  well  . 

With  the  liquid  filled  shell,  however,  there  is  no  a  priori  assurance 
that  an  epicyclic  fit  could  be  made.  This  is  so  for  several  reasons. 
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For  the  solution  of  the  differential  equation  of  the  yawing  motion, 
applicable  to  the  data  in  the  spark  range,  it  is  assumed  that  over  the 
relatively  short  time  of  flight,  all  the  aerodyanmic  coefficients  are 
constants.  The  only  dependent  variable  which  is  permitted  to  vary  slowly 
is  the  nondimens ional  axial  spin.  Its  slow  variation  is  taken  care  of 
by  a  perturbation  technique.  However,  for  liquid  filled  shell,  the  axial 
spin  varies  much  more  rapidly,  at  least  initially,  than  is  the  case  for 
a  rigid  shell.  Moreover,  the  liquid  produces  additional  forces  and 
moments  which  may  not  be  constants.  Both  of  these  effects  may  change 
the  epicyclic  character  of  the  yawing  motion. 

To  test  the  effect  of  the  observed  variation  of  axial  spin  on  the 
yawing  motion,  analogue  computer  runs  were  made  for  two  liquids:  water 
And  mercury.  In  spite  of  the  considerably  faster  initial  decay  of  the 
axial  spin  with  water,  than  would  be  the  case  if  the  filler  were  rigid, 
no  significant  effect  on  the  character  of  the  epicycle  could  be  detected. 
With  mercury,  the  story  is  different.  Because  of  the  extremely  rapid 
initial  decay  of  the  axial  spin  in  this  case,  the  yawing  motion  becomes 
divergent  but  finally  settles  down  to  a  normal  damped  motion  by  mid¬ 
range.  It  should  be  kept  in  mind  that  these  investigations  pertain 
only  to  the  effect  on  the  yawing  motion  of  postulated  variable  axial 
spin  and  have  nothing  to  do  with  possible  other  effects  the  liquid  may 
or  may  not  have  on  such  a  motion. 

It  was  pleasantly  surprising,  therefore,  that  for  all  liquids  provided 
that  the  yaw  was  not  excessive,  i. e.  B  —  10°,  the  observed  yawing  motion 
could  be  represented  quite  well  by  the  usual  epicycle  of  the  formal 

\  =  K1q  exp  (-0^  +  ijtf'-^P  +  k2Q  exp  (-C^  +  i02*)p 
where  \  complex  yaw 

2  are  the  nutational  and  precessional  yaw  damping  rates  respec¬ 
tively 

0n  £  axe  the  circular  rates 

p  the  distance  along  the  trajectory  in  calibers. 
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For  a  rigid,  shell,  the  yaw  damping  rates  are  given  by  the  following 
expressions: 
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The  accuracy  of  the  epicyclic  representation  of  the  observational  data 
varies  somewhat  depending  on  the  magnitude  of  the  yaw  and  the  nature  of 
the  liquid.  For  water,  for  example,  the  least  square  fit  of  an  epicycle 
to  the  observational  data  is  characterized  by  the  standard  deviation  of 
about  7  minutes  of  arc.  This  is  comparable  with  the  accuracy  of  the 
data  and,  hence,  can  be  considered  as  satisfactory.  A  more  visual  sense 
of  the  goodness  of  fit  can  be  obtained  from  the  plots  on  semi-log  paper 

A  4-"U  n  n  1  n  n  nts  A  A  -P  1  rr\T  T  n  ««  -!v-l/-l4-  -?  n4  n  A  rt  rt  /-l  VI/3  -|  ■*->  t~+ 

ui  biic  euujjx  x  uuuetj  ui  eauxx  uum^uncub  Ui  yaw  a&aii uid  tautc  •  Av-^uiuiug, 

to  the  linearized  theory,  the  logarithmic  decrement  of  the  amplitude  of 
yaw  should  be  a  constant,  i.e.,  on  semi-log  paper  the  plot  of  the  amplitude 

iro  o4onoo  oViauI  r}  ho  ^cno-pocont  ohl  o  Vrw  a  a  -4  v*q  -f  rrVi  +■  1  ^  ru=»  .Q<=»tr<=»r»o  1  c:i  i  r»Vi  'nl 
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are  reproduced  in  Figure  2.  The  plotted  points  are  from  a  preliminary  hand 
reduction  prior  to  the  machine  reduction  and  the  application  of  differential 
corrections.  On  the  whole,  straight  lines  do  fit  the  data  reasonably  well. 
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One  can  now  go  farther  and.  complete  the  analysis  of  the  yawing  notion 
by  computing  what  one  may  call  "effective  aerodynamic  coefficients",  i.e., 
from  the  observed,  yaw  damping  rates  compute  H,  T,  and  D.  It  is  to  be 
understood,  of  course,  that  for  a  liquid  filled  shell,  these  are  not  the 
true  aerodynamic  coefficients  but  are  simply  certain  parameters  of  the 
epicycle.  Only  if  it  can  be  demonstrated  that  their  numerical  values 
are  the  same  as  for  an  equivalent  rigid  shell  can  they  be  called  aerodynamic 
coefficients.  Then  the  inference  can  be  drawn  that  the  liquid  does  not 
contribute  measurable  amounts  to  corresponding  aerodyanmic  forces  and 
moments. 
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5.  THE  EFFECT  OF  WATER 

The  water,  as  a  filler,  was  tested  most  extensively.  The 
"effective  aerodynamic  coefficients"  and  the  circular  frequencies, 
for  this  liquid,  are  given  in  Figures  3  and-  4.  The  scatter  in 
the  lift  coefficient  is  to  be  expected  because  it  is  determined  from 
the  swerving  motion  which  method  is  less  accurate  than  that  using 
two  center  of  mass  positions.  The  scatter  in  the  other  coefficients 
is  somewhat  larger  than  one  would  expect  if  the  shell  was  rigid.  Thi6 
might  reflect  the  fact  that  the  fitted  epicycle  is  a  fairly  coarse 
approximation  to  the  actual  yawing  motion.  Also,  there  appears  to  be 
a  systematic  difference  in  the  observed  "coefficients"  between  lower 
and  higher  percentages  of  fill.  To  test  the  significance  of  this 
difference,  the  data  were  divided  into  three  groups.  The  first  group 
consisted  of  all  rounds  from  0-35  percent  of  fill;  the  second  from 
50  to  90  percent  of  fill  and  the  third  group  consisted  of  rounds  of 
100  percent  fill. 

The  following  table  gives  the  average  values  for  each  group  to¬ 
gether  with  its  standard  deviation  and  the  significance  of  the  differ¬ 
ence  of  100  percent  and  50  -  90  percent  groups  relative  to  the  0-35 
percent  group. 

For  this  particular  shell  and  this  particular  liquid,  the  principal 
effect  of  the  liquid  filler  seems  to  be  a  decrease  in  the  effectiveness 

of  the  yaw  damping  moment,  C  +  C  .  Other  effects  are  relatively 

.  q  a 

minor . 

The  yaw  damping  rates  are  the  logarithmic  decrements  of  the 
amplitudes  of  the  precessional  and  nutational  arms  of  the  epicycle. 

In  ballistic  practice,  these  are  defined  in  such  a  way  that  a  positive 
yaw  damping  rate  corresponds  to  shrinking  of  the  amplitude  and  a  nega¬ 
tive  rate  to  its  growth.  Thus,  a  dynamically  stable  shell  should  have 
a  positive  yaw  damping  rate  for  each  component  of  yaw. 

Figure  5  shows  the  yaw  damping  rate  for  water  as  the  liquid  filler 
for  various  percentages  of  fill.  As  was  expected,  the  liquid  affects 
only  the  nutational  yaw  damping  rate.  The  precessional  component  is 
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TABLE  I 


Liquid.  (H^O)  Filled  Shell 
Average  "Aerodyarimic  Coefficients" 


Coefficients 

Group  I 

0  -  55$ 

Group  III 

100$ 

Significance 

Group  II 

50  -  90$ 

Significance 

.0428  +  . 0009 

.0451  +  .0011 

n.  s.* 

.0417  +  .0014 

Signif. 
on  5$ 
level 

.0424  +  .0016 

.0459  +  -0024 

n.  s. 

.0428  +  .0025 

n.  s . 

.0057  +  *0010 

.0029  +  .0019 

n.  s. 

.0015  +  .0021 

s  • 

C,„  +  CM 

Mq  Ma 

-.no  +  .on 

: 

-.079  +  .008 

V.S.* 

+.029  +  .046 

| 

v.s. 

2.506  +  041 

2.502  +  .047 

EL.  s. 

2.260  +  .052 

Signif. 
on  2$ 

% 

.187  +  .006 

. 192  +  . 004 

EL.  S. 

.191  +  .006 

:n.s. 

*  n.  s.  =  not  significant 
v.s.  =  very  significant 


practically  unaffected.  What  was  unexpected  is  a  broad  zone  of  fill 
conditions,  from  about  40  percent  to  90  percent,  for  which  this  com¬ 
ponent  of  yaw  is  divergent.  The  shell,  therefore,  for  these  conditions 
of  fill,  is  dynamically  unstable.  There  is  a  curious  scatter  around 
40  -  50  percent  of  fill  which  remains  unexplained. 

It  is  worthy  of  note  that  Stewartson's  theory  predicts  instabil¬ 
ity,  for  this  case,  in  a  narrow  zone  of  4-3  percent  to  48  percent 
of  fill.  This  is  the  only  zone  of  instability  predicted.  The  fact 
that  at  50  percent  of  fill  the  divergence  of  yaw  is  very  mild  suggests 
that  the  fluid  may  be  approaching  a  condition  of  spin  at  which  Stewartson»s 
theory  may  begin  to  be  applicable.  With  full  spin  of  the  liquid  at  50 
percent  fill  the  shell  should  be  stable.  Since  our  fill  conditions 
were  not  precisely  controlled,  in  that  they  were  predicated  on  the  average 
size  cavity,  it  might  well  be  that  the  four  discordant  points  actually 
lie  in  the  above  zone,  45  -  48  percent,  although  nominally  were  designated 
as  40  percent  and  50  percent  respectively.  If  this  is  so,  then  the  dis¬ 
cordant  points,  showing  larger  divergence  of  yaw  than  those  at  50  percent 
fill,  might  be  considered  as  a  suggestive  confirmation  of  theoretical 
prediction. 
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4.  THE  EFFECT  OF  OTHER  LIQUIDS 

In  addition  to  water,  several  other  liquids  were  used  as  fillers 
in  a  limited  firing  program.  These  liquids  were: 


Liquid 

Sp.  gr. 

V  c.  s. 

Mercury  (Hg) 

15.5 

0.1 

Tetrabromoe thane  (CHBr^) 

5-0 

5.4 

Glycerin 

1.26 

1000 

The  yaw  damping  rates  for  these  liquids  are  given  in  Figure  6. 

Average  values  for  water  are  included  for  comparison. 

Several  interesting  features  are  to  be  noted  in  this  figure.  The 
severity  of  instability  increases  with  an  increase  in  the  specific  gravity 
of  the  liquid.  Thus,  at  "JO  percent  of  fill,  where  the  growth  of  the 
nutational  component  of  yaw  appears  to  be  most  severe,  for  water,  this 
component  of  yaw  will  double  in  about  5460  calibers  of  travel  (226  feet), 
for  CHBr^  in  1150  calibers,  and  for  mercury,  in  about  550  calibers.  The 
pattern  of  instability  versus  percent  of  fill  remains,  however,  essentially 
similar  as  if  the  water  curve  was  pulled  down  roughly  in  proportion  to 
the  ratio  of  the  specific  gravities  of  the  liquids  concerned. 

The  analysis  of  mercury  filled  rounds,  because  of  very  rapid 

divergence  of  yaw,  is  less  certain.  That  is  why  the  symbols  in  Figure  6 

are  placed  in  parentheses.  It  is  to  be  recalled  that  analogue  computer 

run  for  a  completely  mercury  filled  cavity  showed  an  initial  divergence 

of  yaw  due  to  very  rapid  spin  decay.  This  divergence,  however,  is  very 

mild  in  comparison  with  the  observed  divergence  for  a  100  percent  full 

-4 

case.  The  former  rate  of  growth  of  yaw  is  less  than  0.1  x  10  per 
caliber  of  travel;  the  observed  divergence  is  10  x  10  or  one  hundred 
times  greater.  Therefore,  within  the  accuracy  of  the  analysis  of  these 
rounds,  the  effect  of  the  spin  decay  on  the  yawing  motion  of  the  shell 
has  been  neglected. 

The  similarity  between  the  three  curves  may  be  somewhat  fortuitous, 
occurring  only  at  very  early  stages  of  flight  when  the  spins  of  these 
three  liquids  are  far  removed  from  the  full  spin  condition.  Nonspinning 
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liquid  produces  no  detrimental  effect  on  the  dynamics  of  the  shell. 
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produces  no  detrimental  effect  as  predicted  by  Stewart son  and  verified 
by  our  firings  with  glycerin.  Therefore,  instability  for  this  loading 


condition  commences  s,t  some  condition  of*  fluid  motion^  t)6  it  spin  or 

secondary  flows  or  both,  and  gradually  ceases  at  a  spin  distribution 

less  than  full,  because  theoretically  it  requires  a  relatively  long 

time  for  the  liquid  to  reach  full  spin.  The  yaw  damping  curves  for  the 

three  liquids,  H^O,  CHBr^,  and  Hg,  Figure  6,  correspond,  therefore, 

because  of  their  different  viscosities,  to  sonewhat  different  states  of 

fluid  motion.  At  later  times,  the  observed  relationship  between  their 

/  — 

yaw  damping  rates  will  be  altered. 


Glvcerin.  because  of  its  verv  high  viscositv.  attains  rigid  bodv 
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considerable  fraction  of  the  full  spin  is  attained  while  the  shell  is  in 
the  gun.  Hence,  Stewartson's  theory,  for  an  inviscid  fluid,  may  con- 

**  /•*  -t  ir  n  ymti!  tr  +/\  4-Vi  ■?  n  Art  o  rk  A  rt  -r\  vkrti  r-i  Ai  ini  ir  nuAm  4-  -I  Art  0-4-  atta  v»+-  nrtw  -rt-vlrt 
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diets  instabilities  only  between  45  percent  and  48  percent  of  fill  and 
nowhere  else.  The  shell  have  been  tested  at  90  percent,  70  percent 
and  48  percent  fill 

percent.  This,  therefore,  may  be  considered  as  a  verification  of 
Stewartson's  theory. 


As  Figure  6  shows,  the  shell  is  unstable  only  at  48 


5.  AXIAL  SPIN 


As  previously  mentioned,  dissimilar  pins  in  the  base  of  the  shell 
permit  determination  of  its  axial  spin  history  in  the  range.  Some  typical 
spin  histories  for  an  empty  shell,  100  percent  filled  with  water,  glycol, 
thetrabromoethane,  and  mercury  are  shown  in  Figure  J.  It  is  seen  that 
the  spin  decay  of  an  empty  shell  is  very  modest;  it  decreases  in  the  first 
120  milliseconds  of  flight  by  about  1.6  percent.  If  the  cavity  was  filled 
with  a  rigid  substance,  inertially  equivalent  to  water,  the  spin  would 
decrease  only  1.5  percent.  However,  with  water  it  decreases  by  5*6  percent, 
with  thetrabromoethane  by  15  percent,  and  with  mercury,  in  the  first  60 
milliseconds  of  travel,  by  25  percent.  However,  even  with  mercury  the  shell 
remained  gyroscopically  stable  in  spite  of  a  substantial  loss  of  its  axial 
angular  momentum. 

The  relatively  rapid  spin  decay  of  the  liquid  filled  shell  is  due,  of 

course,  to  a  continuous  transfer  of  shell  angular  momentum  to  the  liquid. 

The  two  are  coupled  by  the  friction  at  the  walls  of  the  cavity  and  the  radial 

velocity  gradient,  thus  established,  transmits  angular  momentum  to  the.  liquid. 

The  efficiency  of  this  process  depends  on  whether  the  momentum  is  transferred 

by  laminar  o t  the  turbulent  shearing  stresses  at  the  wall.  The  latter  is 

many  times  more  efficient  than  the  former.  However,  initially,  the  boundary 

layer  is  probably  laminar  because  of  the  relatively  low  Reynolds  number. 

If  one  defines  the  Reynolds  number  as 

(V  -  V)B 

O 


where  VQ  is  peripheral  velocity  of  the  shell 

V  is  the  velocity  at  the  edge  of  the  boundary  layer  (V  =  0  initially) 
and  6  is  the  thickness  of  the  boundary  layer,  then  initially  Re  is  small 
because  of  small  8.  It  increases  to  its  maximum  value 


Re 


V  a 
o 


max 


and  then  decreases  because  of  the  decrease  in  the  difference  between  a>c 
and  oo  as  the  liquid  approaches  rigid  body  fotation;  i.e., 
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Re 


For  the  condition  of  our  experiments  with  water,  the  maximum  Re  is 
about  one -half  million.  This  value  is  reached  in  about  two  seconds  of 
time  of  flight.  In  the  range,  therefore,  where  observations  extend 

only  to  0.1  seconds,  the  maximum  Reynolds  number  is  of  the  order  of  30*000. 

12 

Schlichting  has  studied  the  stability  of  the  laminar  boundary  layer 
inside  a  spinning  cylinder  whose  spin  is  started  inqpulsively  from  rest,  a 
situation  not  dissimilar  to  ours.  He  has  shown  that  subject  to  small  dis¬ 
turbances,  the  laminar  boundary  layer  is  stable  up  to  Reynolds  number 
of  66,000,  which  is  the  minimum  critical  Reynolds  number.  Therefore,  in  our 
case,  the  transition  to  turbulence  is  not  likely  to  occur  during  the 
observed  flight  of  0.1  seconds. 

In  view  of  these  considerations,  it  was  assumed  that  the  fluid 
motion  can  be  represented  by  a  two  dimensional  laminar  flow.  A  theoretical 
analysis  of  this  problem  was  made  by  Scott,  reference  11.  Figure  8  shows 
that  for  water,  the  two  dimensional  laminar  flow  hypothesis  is  inadequate. 

The  observed  spin  decay  can  be  represented  as  a  product  of  two  expo¬ 
nent  iais 

-Px  -f 
o>  =  0)  e  e 
o 

where  the  first  exponential  is  due  to  the  air  torque  and  the  second  due 
to  the  liquid. 

Using  the  usual  laminar  boundary  layer  approximations,  it  can  be 
shown  that  laminar  shear  stress  on  the  walls  of  the  rotating  cylinder, 
started  impulsively  from  rest,  is 


where  a)  is  the  spin  of  the  cylinder,  a  is  radius  of  the  cavity  and  p,  v 
are  the  density  of  the  liquid  and  its  kinematic  viscosity  respectively. 
The  torque  due  to  the  liquid,  therefore,  is 

T.  =  a  .  tit  .  S 
i  W 
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where  S  is  wetted  area.  The  differential  equation  governing  the  spin  decay 

_  _a  xi _ _■» _ i  i _ _  t _ _ „ 

Oi  brie  snexx  utiu  ue  wriiucu  etc> 


A  a>  =  -  T  ,  -  T,  . 

air  xiquia 

2 

=  -  Pa  S  U  cq  C  T, 

P  " 


where  A  is  the  polar  moment  of  inertia  of  empty  shell.  From  this,  it 
follows  that 


f  = 


r 


j 


.  T< 
=  1  jc 


dt 


and 


p 


d2.qf! 


~ZK 


p 


For  the  laminar  boundary  layer,  using  the  expression  for  the  liquid 
torque  and  performing  the  integration,  the  exponent  "f"  can  be  written  as 


f  =  cx  l/T 


where  is  a  computable  constant; . 


iiiio  CApiCDDiUU 


of  Figure  8 ,  which  was  computed  by  a  different  method  in  reference  ll,  very 
well.  This  is  not  surprising.  The  thickness  of  the  laminar  boundary 


layer,  at  the  end  of  the  observed  flight,  is  only  about  10  percent  of  the 
radius  of  the  cavity.  Therefore,  the  boundary  layer  solution  is  a 
reasonable  approximation  to  the  more  exact  solution.  At  later  times,  the 
two  solutions  should  become  progressively  divergent. 


For  the  turbulent  boundary  layer,  the  unsteady  solution  is  unknown. 
However,  in  order  to  get  some  idea  of  what  the  spin  decay  would  be  if  the 
boundary  layer  was  turbulent.  Dr.  Sternberg  suggested  the  use  of  the 
momentum  integral  with  an  assumed  velocity  distribution.  The  turbulent 
shear  stress  at  the  wall  is 


8 

ar 

Tw  =  “a —  I  p  u  a  y 

''o 

if  one  assumes  the  velocity  distribution 


iL  =  n  _  l) 

U  6 

o 


the  shear  stress  becomes 


rW  =  B  p  Uo  dt 


The  toraue  T .  =  a.  S  .  txt 
■*  l  W 

and  the  "f"  value 


r  T.  s 

f’J  ^dt=J 

o  o 


s*  aSx 


aS  i 
o 


rttl  f\(+\ 

Zo  (5) 

Aob 


if  ob  is  assumed  constant. 


To  get  6(t)  we  proceed  as  follows: 


1  _  rr  2  „  I 


a  K  “o  at 


Using  flat  plate  results: 


and  eliminating  x  between  eqs  (7)  and  (8)  one  obtains  Cp  =  g(6). 
Substituting  in  eq  (6)  end  integrating  one  gets 


6(t)  »  .0639  (Uq) 


1/5  4/5  4/5  4/5 


24 


Finally,  substituting  (9)  into  (5)  one  gets  the  desired  result.  The 
turbulent  case  of  Figure  8  was  computed  in  this  manner.  The  accuracy 
of  these  approximations  is  unknown  but  the  results  should  be  adequate 
for  illustrative  purposes.  It  may  be  of  interest  to  note  that  if  one 


used  the  flat  plate  local  turbulent  skin  friction  coefficient  directly. 


T  =  C 
W  f 


a  = 


.05' 


BK 


and  substituting  for  x  =  U  t  in  R ,  one  obtains  a  similar  expression 
w  o  “ 

for  "f"  of  the  form 


Only  the  numerical  factor  is  1.55  times  smaller  than  the  numerical  factor 
obtained  by  the  momentum  method.  Figure  8  shows  that  neither  the  laminar 
nor  turbulent  boundary  layers  produce  the  type  of  torque  which  is  necessary 


to  account  for  the  observed  spin  decay.  The  observed  curve  lies  in  be¬ 
tween  these  two  extremes.  Nevertheless,  the  observed  spin  decay  can 


better  be  represented  by  the  vt 
layer  is  laminar. 


law  suggesting,  thereby,  that  the  boundary 


Scott  has  found  that  it  is  possible  to  represent  the  observed  spin 
decays  of  various  liquids,  with  the  two  dimensional  laminar  hypothesis, 
by  using  fictitious  viscosities.  Thus  for  water,  an  "effective"  viscosity 
52  times  greater  than  its  natural  viscosity,  gives  an  excellent  representa¬ 
tion  of  the  Observed  spin  decay;  for  glycol,  however,  its  natural  viscosity 
suffices,  Figure  9.  The  following  fictitious  viscosities  represent  the 
observed  spin  decays  of  various  liquids  * 


Liquid 

v  c.s. 

0 

v/v( 

Hg 

0.11 

290 

H2° 

1. 

52 

CHBr2 

3.4 

14 

Glycol 

660 

i 

If  the  boundary  layer  is  not  turbulent,  and  it  does  not  appear  to 

be,  then  the  "effective"  high  viscosities  must  be  accounted  for  by  some 

mechanism  which  increases  the  laminar  shear  stress  at  the  wall.  Such 

mechanism  may  be  provided  by  the  secondary  flows  in  the  cavity.  It 

is  well  known  that  in  an  inclosed  cylinder  filled  with  a  liquid  and 

started  to  spin  impulsively,  secondary  flow  develops  as  a  result  of  end 
13 

effects  .  Such  circulation  may  well  steepen  the  velocity  gradient 
at  the  wall  and,  hence,  effectively  increase  the  laminar  shear  stress. 

In  order  to  examine  this  hypothesis  in  greater  detail,  an  "impulsive 

Ik 

spin  generator"  was  designed  and  built'*'".  It  consists  of  a  transparent 
cylinder,  geometrically  homologous  to  the  shell  cavity,  mounted  on  a 
rotating  table  which  could  be  instantaneously  engaged  to  a  rotating  fly 
wheel  operated  by  an  air  turbine.  Two  magnetic  pickup  coils  measure 
the  speed  of  rotation.  If  spun  at  7000  rpm,  the  apparatus  correctly 
simulates  the  ratio  of  centrifugal  to  deceleration  forces  of  the  shell 
in  free  flight. 
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of  a  turbulent  boundary  layer .  Tufts  were  unsuccessful.  The  use  of 
"aluminum  lining11  powder,  flaky  bits  of  aluminum  which,  in  suspension, 
align  themselves  in  the  direction  of  shear  surfaces,  is  a  well-known 
technique  for  visual  detection  of  turbulence.  The  results  were  negative; 
the  flow  appeared  to  be  laminar. 

To  detect  and  measure  the  secondary  flows,  a  method  was  devised 
employing  a  homogeneous  suspension  of  particles  in  the  fluid.  The 
particles  were  black  ethvl  cellulose  (0  =  l.ll)  suspended  in  an  aqueous 

A -  W  V  w  —  — 

solution  of  gyleerin,  the  solution  having  a  kinematic  viscosity  of  .062 
Stokes.  Their  motion  was  photographed  with  a  high-speed  camera.  Fuller 
account  of  these  experiments  will  be  given  elsewhere.  Here,  only  a  few 
pertinent  results  will  be  mentioned. 


When  the  film  was  projected,  the  existence  of  secondary  circulation 


was  at  once  apparent • 
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upper  and  lower  cells.  In  each  cell,  a  doughnut-like  vortex  was  generated, 
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carrying  fluid  particles  from  the  periphery  toward  the  axis  and  up. 
in  the  upper  part  of  the  cylinder,  and  down,  in  the  lower  part.  Figure 
10  shows  the  trajectories  of  some  of  the  particles.  These  show  clearly 
the  two-cell  structure  of  the  fluid  notion  and  a  strong  indication  of  a 
development  of  a  doughnut -like  vortex  in  each  cell. 

Distribution  of  vorticity  throughout  the  fluid,  at  two  instances  of 
time,  is  shown  in  Figure  11.  Corresponding  theoretical  distribution  for 
laminar  flow,  no  end  effects,  and  without  the  secondary  flows,  is  given 
for  comparison. 

In  the  solution  of  the  two-dimensional  laminar  flow  problem^, 
viscosity  and  time  appear  only  as  a  product,  vt.  Therefore,  the  observed 
distributions  of  vorticity,  at  various  peripheral  speeds,  can  be  roughly 
approximated  by  the  theoretical  distributions  either  with  fictitious 
viscosities  and  real  times  or  natural  viscosity  and  later  times.  The 
following  tables  give  the  ratios  of  the  fictitious  to  natural  viscosity 
for  various  cases.  With  the  fictitious  viscosity,  theoretical  curves 
could  be  approximately  fitted  to  the  observed  points. 

TABLE  II 

Matching  the  observed  distributions  of  vorticity  by  fictitious 
viscosities  and  real  times. 

v/vQ 


Observed  Time,  t  sec. 

oe 

II 

-0 

H 

1434 

3065  rpm 

1 

l 

1.5 

2 

2 

1.5 

2 

5 

5 

1-7 

2.3 

4 

4 

2.5 

The  results  show  that  the  greater  is  the  peripheral  speed,  the 
higher  is  the  "effective"  viscosity  and,  hence,  the  effective  shear 
stress  at  the  wall.  However,  the  details  of  the  mechanics  whereby  the 
secondary  circulation,  of  the  observed  nature,  increases  the  laminar 
shear  at  the  wall,  are  not  clear. 
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If  the  trend  of  the  increase  in  the  effective  viscosity  with  the 

increase  in  the  peripheral  speed,  as  shown  by  table  experiments, 

continues,  then  one  would  expect  considerably  higher  ratios  of  v/v 

5  ° 

for  the  20mm  shell  whose  spin  is  about  10  rpm.  The  following  empirical 
relationship  seems  to  tie  the  20mm  shell  and  the  rotating  cylinder 
results  in  a  reasonable  way,  although  its  physical  significance,  if  any, 
is  not  at  all  clear.  The  ratio  of  the  effective  viscosity  to  the  natural 
viscosity  for  various  liquids  and  spins  can  be  represented  by  the  follow¬ 
ing  relationship: 


A 
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where  viscosity  of  the  liquid  in  centi-stokes,  and  v  and  vr  are  the 
peripheral  velocities  of  the  rotating  cylinder  (table  experiments)  and 
the  20mm  she] I  respectively.  A  comparison  between  the  "observed."  values 
and  the  empirical  formula  above  is  shown  in  the  following  table: 


cvyumt  kJn^-L-L. 


Liquid 

V  cm/sec. 

F  -Z 

v  c.  s. 

0 

v/vq  obs. 

v/vq  emp. 

Hg 

7.6  x  ioJ 

0.11 

290 

290 

w  o 

2 

CHBr2 

7=6  x  105 
_  ^  .  ~3 

(.OX  XU' 

1.0 

5*4 

52 

1  \. 

-L*f 

32 

t  r\ 

XVJ 

Glycol 

7.6  X  103 

660 

1 

1 

Rotat ing 

v  cm/sec. 

V  c.  s. 

n 

v/ v  obs . 

•  n 

v/v  emp. 

Cylinder 

O 

H20  +  Glycerin  603.4 

6.2 

3.0 

2.6 

H20  +  Glycerin  282.3 

6.2 

1.9 

2.0 

KgO  +  Glycerin  92.7 

6.2 

1.4 

1.4 

Tne  agreement  is  not  bad  and  the  empirical  formula  can,  perhaps,  be 
used  to  estimate  an  early  history  of  the  shell fs  axial  spin  decay. 

However,  even  with  the  enhanced,  fictitious  viscosity,  it  takes  a 
relatively  long  time  for  a  low  viscosity  fluid,  such  as  water,  to  acquire 
full  spin  and  reach  a  state  at  which  Stewartson!s  theory  becomes  applicable. 
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Meanwhile,  if  the  shell  is  dynamically  unstable  during  the  transition  to 
the  full  spin  regime,  the  yaw  may  become  prohibitively  large  rendering  the 
shell  useless  in  practice.  To  test  the  effective  duration  of  the  transient 
regime  and  its  effect  on  the  theoretically  predicted  stable  flight  (when 
the  liquid  is  fully  spinning),  the  20mm  shell  were  fired  at  longer  ranges. 
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6.  LONG  RANGE  FIRINGS 


Six  groups  of  20mm  shell  were  prepared,  five  shells  in  each  group, 
for  firing  at  a  vertical  target  located  600  yards  from  the  muzzle.  The 
shell,  identical  with  those  fired  in  the  spark  range,  were  loaded  with 
various  percentages  of  fill  with  two  types  of  liquids:  water  and  the- 
trabromoethane,  in  accordance  with  the  following  schedule. 


H20  (Sp.  Gr.  =  1) 
Groups  Fill,  percent 

1  100 

2  70 

5  50 


Groups 

4 

5 

6 


CHBr  (Sp.  Gr.  =  3.0) 
rill,  percent 

90 

70 

50 


The  shells  were  to  be  fired  from  a  standard  gun  (twist  of  rifling 
1:25)  at  a  muzzle  velocity  of  2J00  fps.  However,  due  to  an  erroneous 
loading,  the  muzzle  velocity  turned  out  to  be  closer  to  2600  fps  than  to 
27OO  fps.  This  small  difference  in  the  muzzle  velocity  is  unimportant 
for  the  purposes  at  hand. 

The  experiment  consisted  of  measuring  the  velocity  by  two  velocity 
screens  42  apart,  with  the  first  screen  located  26*  from  the  muzzle,  and 
the  usual  time  of  flight  screen  at  600  yards.  A  paper  screen  on  the 
second  velocity  screen  gave  some  indication  of  the  initial  yaw.  All 
initial  yaws  appeared  to  be  small. 


The  analysis  of  the  firings  consisted  in  the  evaluation  of  an 
"effective"  ballistic  coefficient  over  about  one  second  of  the  time  of 
flight. 
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The  results  are  given  in  the  following  table: 

TABLE  III 

Ballistic  coefficients  of  water  filled  20mm  M56  shell  obtained  from 


4-  A  -PI  4  rrV.4- 

OX1UC  U1  IIJ.51IU 

firings  to  600  yards. 

Percent  of 

Weight  of  Shell 

Fill 

Grains 

C 

C  -  Expected 

i  00 

xvv 

1128 

•  27k 

- ,  . 

.275 

1131 

r\t-»  O 

.  d  f  U 

1129 

ro 

-0 

1129 

.273 

70 

IO87 

.216 

.266 

IO89 

.218 

IO87 

ORO 

. 

1090 

ro 

0 

pzr\ 

1088 

.266 

=  265 

1100 

.256 

1080 

.259 

1077 

- II 

.260 

1092 

.238 

The  average  cavity  requires  7=8  milliliters  to  fill,  corresponding 
to  a  weight  of  water  of  about  120  grains.  Apparent  Inconsistencies  in 
the  above  weights  arose  from  variations  in  the  weights  of  the  empty  shell. 

What  was  the  expected  behavior  of  these  shell  as  inferred  from  the  spark 
range  results,  illustrated  in  Figure  5?  The  100  percent  filled  shell  are 
stable  and,  hence,  should  fly  well.  The  consistency  of  the  ballistic 
coefficients  and  their  agreement  with  the  expected  value  verifies  this. 
According  to  Figure.  5  the  shell,  in  the  transient  regime,  with  70  percent 
water  filled  cavity,  is  unstable.  The  degree  of  instability  is  such 
that  the  yaw  doubles  every  5460  calibers  of  travel.  Since  the  distance 
to  the  600  yard  target  is  about  28,000  calibers  if  this  divergence  of  yaw 
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persisted  throughout  the  flight,  the  yaw  on  the  target  would  he  about 
times  as  great  as  the  initial  yaw.  Even  if  the  initial  yaw  is  of  the 
order  of  one  degree  or  so,  the  final  yaw  would  be  prohibitively  large. 

However,  if  the  "full"  spin  is  reached  at,  say,  one-third  of  the  range, 
then  the  shell  will  become  stable,  as  was  demonstrated  by  our  firings  with 
glycerin,  and  from  this  point  on  the  yaw  will  decay.  If  this  history 
of  yaw  is  followed,  on  would  expect  that  the  average  ballistic  coefficient 
will  be  decreased  by  about  1.6  percent  where  6q,  the  initial  yaw,  is 
in  degrees.  The  50  percent  filled  shell  were  also  expected  to  show  lower 
ballistic  coefficient.  The  divergence  of  yaw,  at  this  fill,  is  milder 
than  at  70  percent  and  the  liquid  should  reach  "full"  rotation  sooner. 

However,  the  behavior  of  shell  at  this  percentage  of  fill  is  on  the 
ragged  edge  because  at  48  percent  fill,  the  shell  is  predicted  to  be  un¬ 
stable  as,  again,  was  verified  by  the  spark  range  firings  with  the  glycerin. 
Therefore,  in  a  group  of  shell  fired  with  this  percent  of  fill,  there  might 
be  an  occasional  "maverick"  due  to  inadequately  controlled  loading  conditions. 
One  of  the  rounds  in  the  50  percent  fill  group  appears  to  be  in  this  category. 

The  story  of  firings  with  thetrabromoethane  is  short.  Not  a  single 
shell  of  fifteen  fired  hit  the  target  in  spite  of  the  valiant  efforts  to 
do  so.  All  fell  short  even  when  the  gun  was  elevated  to  shoot  40  feet 
above  the  target. 

The  rate  of  divergence  of  the  nutational  component  of  yaw  for  CHBr^ 
is  about  three  times  as  large  as  for  water.  Therefore,  at  corresponding 
times  in  flight,  the  yaw  levels  for  these  rounds  must  have  been  considerably 
larger.  This  demonstrates,  rather  dramatically,  the  importance  of  the  un¬ 
stable  transient  regime  for  heavier  liquids. 
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7*  BAFFLED  ROUNDS 

The  analysis  of  the  axial  spin  of  the  shell  shows  that,  during  the 
observing  period  in  the  firing  range,  all  fluids,  with  the  exception  of 
such  high  viscosity  liquids  as  glycerin,  v  =  1000  c.s.,  or  glycol, 
v  =  660,  continue  to  absorb  axial  angular  momentum  from  the  shell.  Hence, 
they  are  in  an  unsteady  state,  only  gradually  reaching  rigid  body  rotation. 
It  is  clear,  therefore,  that  the  dynamic  instability  of  the  nutational 
component  for  a  broad  range  of  fill  conditions  is  somehow  associated  with 
this  transient  phase  of  liquid  motion.  The  mathematics  of  this  problem 
is  extremely  complicated.  Certain  tentative  approaches  are  being  made 
and  the  solution  may  eventually  be  found. 

However,  from  the  engineering  point  of  view,  it  would  appear  that 
this  transient  unstable  phase  can  be  eliminated  by  forcing  the  liquid 
to  acquire  full  spin  at  the  muzzle.  This,  of  course,  can  be  done  by 
inserting  baffles  into  the  cavity.  It  is  a  crude  approach  but  it  was 
tried  to  test  the  expectations.  The  baffles  consisted  of  two  cross 
plates  extending  from  the  base  of  the  cavity  to  the  base  of  the  fuze.  The 
following  fillers  were  fired  with  baffles  with  the  following  dynamical 
results: 

Filler  Per  Cent  Motion 


CHBr2 

70 

Stable 

50 

Unstable 

Hg 

100 

Stable? 

90 

Stable 

TO 

Stable 

50 

Unstable 

The  yaw  of  the  100  percent  Hg  round  was  too  small  for  a  reliable  analysis 
but  it  appeared  to  be  stable,  at  least  the  yaw  remained  immeasurably 
small  throughout  the  firing  range . 

It  is  interesting  to  note  that  again  the  50  percent  filled 
rounds  were  unstable  as  predicted  by  Stewartson  for  this  heavy  liquid, 
assuming  his  theory  can  be  applied  to  baffled  rounds.  However,  for 
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these  heavy  liquids,  his  theory  also  predicts  instability  at  90  percent. 
With  baffles,  the  90  percent  filled  Hg  round,  however,  was  stable. 

Thus,  a  doubt  is  cast  on  the  applicability  of  Stewartson' s  theory  to 
baffled  rounds,  at  least  to  the  type  of  baffles  used  in  these  experiments. 
Perhaps  at  50  percent  fill,  the  baffles  do  not  interfere  with  the  fluid 
motion  as  much  as  at  90  percent.  Therefore,  Stewartson*  s  theory  may  be 
valid  at  50  percent  fill  and  not  at  90  percent  fill. 
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8.  SOME  ENERGY  CONSIDERATIONS 


The  fluid  in  a  cylindrical  cavity,  rotating  with  the  full  spin  of 
the  shell,  has  an  infinite  number  of  free  periods  of  oscillations.  These 
frequencies  depend  on  the  geometry  of  the  cavity,  spin  of  the  liquid, 
and  the  amount  of  liquid  in  the  cavity  or  the  percent  of  fill.  The  latter 
establishes  the  free  boundary  conditions.  It  has  been  shown  by  Stewartson 
that  instabilities  of  the  shell-liquid  system  occur  whenever  some  frequency 
of  the  liquid  coincides  with  the  natural  frequency  of  the  solid  shell. 

The  rigid  shell  has  two  natural  frequencies:  precession  and  nutation.  For 
the  usual  shell  design,  only  the  nutational  frequency  lies  in  the  range 
of  fluid  frequencies.  Hence,  only  the  nutational  mode  can  get  into 
resonance  with  the  fluid.  In  comparison  with  the  whole  spectrum  of  fluid 
frequencies,  the  nutational  frequency  is  relatively  low.  Hence,  only  the 
lower  nodes  of  fluid  frequenices  are  usually  of  importance.  These  modes, 
however,  contain  most  of  the  energy  of  the  disturbed  fluid. 

At  resonance,  the  two  oscillating  systems,  the  shell  and  the  fluid, 
can  exchange  energies  very  readily.  The  energy  of  the  disturbed  fluid 
motion  is  transferred  to  the  yawing  motion  of  the  shell  through  work 
done  by  the  pressure  moment  due  to  liquid. 

In  the  transition  phase,  while  the  liquid  is  acquiring  the  spin,  the 
situation  is  probably  different.  Mathematical  solution  of  this  problem 
is  prohibitively  difficult.  Nevertheless,  it  appears  plausible  that, 
with  a  continuous  gamut  of  spins  present  in  the  liquid  during  this  regime, 
the  spectrum  of  discrete  frequencies  of  oscillations  is  absent.  Therefore, 
the  observed  instabilities  of  the  shell  are  not  associated  with  resonance 
conditions  as  was  the  case  in  the  steady  state.  Suitable  transverse 
torques  must  still  be  present  to  keep  the  yaw  growing,  but,  in  the  absence 
of  the  mathematical  model,  the  origin  of  these  torques  is  not  clear. 

It  is  interesting  to  note  how  relatively  little  of  the  energy  of  the 
shell-liquid  system  is  contained  in  the  yawing  motion  of  the  shell.  For 
our  20mm  100  percent  water  filled  shell,  twice  the  kinetic  energy  con¬ 
tents  in  the  various  modes  of  the  angular  motions  are,  approximately, 
as  follows.  The  figures  are  for  the  middle  of  the  observed  flight  in 
the  spark  range. 
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Axial  Spins 
Shell 
Liquid 

Yawing  Motion  (per 
sq.  degree  of  yaw) 
Dissipation 


l*-00  joules 
10  joules 

.02  joules 
16  joules 


To  increase  the  yaw,  therefore,  requires  relatively  little  of  the 
fluid  energy.  As  an  example,  to  account  for  the  observed  decrease  in 
the  yaw  damping  rate  for  the  100  percent  water  filled  shell  (see  Figure 
5),  it  is  necessary  that  the  liquid  moment  be  doing  work  at  a  rate  of 
only  about  k  x  10-^  watts.  This  is  very  small.  Therefore,  in  order  to 
predict  such  an  effect  theoretically  would  require  a  very  precise  knowledge 
of  the  fluid  motion  in  the  transition  regime  and  this  is  not  very  promising. 
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9-  SUMMARY 


For  nonspinning  liquid  filled  shell  both  theory  and  experiment"1 
agree  that  a  liquid  filler  produces  no  detrimental  effect  on  the  dynamics 
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filled  cavity  should,  by  acting  as  an  energy  sink,  increase  the  stability 


of  the  shell  by  increasing  its  yaw  damping  rate. 
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Evidence  of  this  has 
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The  situation  is  quite  different  for  spin  stabilized  shell.  Here, 
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detrimental  effect  on  the  shell  dynamics. 

For  a  shell  with  a  cylindrical  cavity  filled  with  a  liquid,  either 
partially  or  completely,  with  the  liquid  rotating  with  the  full  spin  of 

5 

the  shell.  Stewartson  has  shown  how  to  predict  fill  conditions  at  which 
the  shell  is  dynamically  unstable.  The  designer,  therefore,  may  avoid 
such  loadings. 
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was  necessary  for  Stewartson  to  make  the  assumption  that  the  fluid  is 
inviscid  yet  rotating  with  the  full  spin  of  the  shell.  In  reality,  of 
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state.  Experiments  with  the  20mm  shell  have  shown  the  existence  of  a 
broad  zone  of  loading  conditions  at  which  the  shell  is  dynamically  un¬ 
stable.  This  broad  region  of  instability  is  probably  associated  with  the 
unsteady  phase  of  liquid  motion  in  the  process  of  reaching  rigid  body 
rotation.  The  time  duration  of  this  phase  is  shortened  by  high  viscosity. 
For  a  completely  filled  cylindrical  cavity,  it  is  also  shortened  by  the 
development  of  secondary  flows  within  the  cavity  which  markedly  assist 
in  the  diffusion  of  vorticity  throughout  the  liquid,  thus  expediting 
its  approach  to  the  rigid  body  rotation.  For  partially  filled  cavity, 
the  direct  experimental  evidence  on  the  effectiveness  of  secondary 
flows  is  not  available.  Hie  indirect  evidence,  however,  from  the  ob¬ 
served  spin  history  of  the  partially  filled  shell,  and  their  dynamic 
behavior  when  fired  at  longer  ranges,  suggest  that  the  secondary  flows 
continue  to  be  effective,  at  least  for  higher  percentages  of  fill. 
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Experiments  with  the  20mm  shell  also  indicate  that  the  rate  of  diver¬ 
gence  of  the  nutational  component  of  yaw  increases  with  the  increase  in 
the  specific  gravity  of  the  liquid.  Thus,  for  thetrabromoethane, 
sp.  gr.  =  3,  it  is  roughly  three  times  as  rapid  as  for  water,  sp.  gr.  =  1; 
and  for  mercury  it  is  roughly  thirteen  times  as  fast.  Thus,  in  spite  of 
an  effective  assist  from  secondary  flows  in  reaching  the  steady  state, 
when  Stewartson’s  theory  can  be  used  to  predict  shell  dynamics,  the 
transition  time  may  still  be  long  enough  for  prohibitively  large  yaws  to 
develop.  This  has  been  demonstrated  by  long  range  firings  with  the- 
trabromoethane.  With  mercury,  Plate  1  shows  a  similar  result  only  at 
much  shorter  distance  from  the  muzzle. 

Therefore,  it  is  safe  to  infer  that  for  all  liquids,  with  viscosity 
less  than,  say  500  c.s.,  and  with  specific  gravity>l,  it  is  desirable  to 
suppress  the  transition  phase  by  some  sort  of  baffles.  The  effectiveness 
of  baffles  for  heavier  liquids  has  been  demonstrated  for  liquid  loadings 
in  excess  of  50  percent. 

Whether  the  results  obtained  with  the  20mm  shell  are  scalable  to 
larger  shell  remains  to  be  determined.  The  effect  of  different  geometries 
of  the  cavity  also  needs  examination. 
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